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We use nuclear magnetic resonance (NMR), high-resolution x-ray and neutron scattering to study 
structural and magnetic phase transitions in phosphorus-doped BaFe 2 (Asi_ cc P cc ) 2 . Previous trans¬ 
port, NMR, specific heat, and magnetic penetration depth measurements have provided compelling 
evidence for the presence of a quantum critical point (QCP) near optimal superconductivity at 
x — 0.3. However, we show that the tetragonal-to-orthorhombic structural (T s ) and paramagnetic 
to antiferromagnetic (AF, Tn) transitions in BaFe 2 (Asi_ cc P cc )2 are always coupled and approach 
to Tat « T s > T c (« 29 K) for x = 0.29 before vanishing abruptly for x > 0.3. These results 
suggest that AF order in BaFe 2 (Asi_ x P cc )2 disappears in a weakly first order fashion near optimal 
superconductivity, much like the electron-doped iron pnictides with an avoided QCP. 

PACS numbers: 74.70.Xa, 75.30.Gw, 78.70.Nx 


A determination of the structural and magnetic phase 
diagrams in different classes of iron pnictide superconduc¬ 
tors will form the basis from which a microscopic theory 
of superconductivity can be established m- The par¬ 
ent compound of iron pnictide superconductors such as 
BaFe 2 As 2 exhibits a tetragonal-to-orthorhombic struc¬ 
tural transition at temperature T s and then orders an- 
tiferromagnetically below T/v with a collinear antiferro¬ 
magnetic (AF) structure [Fig. 1(a)] 3 , [4]. Upon hole¬ 
doping via partially replacing Ba by K or Na Pi |U, the 
structural and magnetic phase transition temperatures in 
Bai_ ;E A E Fe 2 As 2 (A = K, Na) decreases simultaneously 
with increasing x and form a small pocket of a mag¬ 
netic tetragonal phase with the oaxis aligned moment 
before disappearing abruptly near optimal supercon¬ 
ductivity mm- For electron-doped Ba(Fei_ ;c T c ) 2 As 2 
(T =Co,Ni), transport [12] 13], muon spin relaxation 
(mSR) QD, nuclear magnetic resonance (NMR) [T514T7] . 
x-ray and neutron scattering experiments [T8ll23l have 
revealed that the structural and magnetic phase tran¬ 
sition temperatures decrease and separate with increas¬ 
ing x [T8ti23] . However, instead of a gradual suppression 
to zero temperature near optimal superconductivity as 
expected for a magnetic quantum critical point (QCP) 
Usuis], the AF order for Ba(Fei_ ;c T c ) 2 As 2 near optimal 
superconductivity actually occurs around 30 K (> T c ) 


and forms a short-range incommensurate magnetic phase 
which competes with superconductivity and disappears 
in the weakly first order fashion, thus avoiding the ex¬ 
pected magnetic QCP [20lf23] . 

Although a QCP may be avoided in electron-doped 
Ba(Fei_ ;E X! E ) 2 As 2 due to disorder and impurity scatter¬ 
ing in the FeAs plane induced by Co and Ni substitution, 
phosphorus-doped BaFe 2 (Asi_ ;r P :E )2 provides an alter¬ 
native system to achieve a QCP since substitution of As 
by the isovalent P suppresses the static AF order and 
induces superconductivity without appreciable impurity 
scattering [241427] . Indeed, experimental evidence for the 
presence of a QCP at x = 0.3 in BaFe 2 (Asi_ ;E P ;E )2 has 
been mounting, including the linear temperature depen¬ 
dence of the resistivity [28], an increase in the effective 
electron mass seen from the de Haas-van Alphen [26] . 
magnetic penetration depth 122113a, heat capacity m, 
and normal state transport measurements in samples 
where superconductivity has been suppressed by a mag¬ 
netic field [32] . Although these results, as well as NMR 
measurements [33], indicate a QCP originating from the 
suppression of the static AF order near x = 0.3, recent 
neutron powder diffraction experiments directly measur¬ 
ing T s and T/v in BaFe 2 (Asi_ ;r P :c )2 as a function of x 
suggest that structural quantum criticality cannot exist 
at compositions higher than x = 0.28 |34|. Furthermore, 
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FIG. 1: (a) The AF ordered phase of BaFe 2 (Asi_ x P cc ) 2 , 

where the magnetic Bragg peaks occur at Qaf = (1,0, L) 
(L — 1,3, •••) positions, (b) Temperature dependence of 
the resistance for the x — 0.31 sample, where RRR — 
T(300 K)/T(0 K) ~ 17. In previous work on similar P- 
doped samples, RRR ~ 13 [28|. (c) The phase diagram of 
BaFe 2 (Asi_ cc P cc ) 2 , where the Ort, Tet, and SC are orthorhom¬ 
bic, tetragonal, and superconductivity phases, respectively. 
The inset shows the expanded view of the P-concentration de¬ 
pendence of T s , T/v and, T c near optimal superconductivity. 
The color bar represents the temperature and doping depen¬ 
dence of the normalized magnetic Bragg peak intensity. The 
dashed region indicates the mesoscopic coexisting AF and SC 
phases. 


the structural and magnetic phase transitions at all stud¬ 
ied P-doping levels are first order and occur simultane¬ 
ously within the sensitivity of the measurements (~0.5 
K), thus casting doubt on the presence of a QCP 34] . 
While these results are interesting, they were carried out 
on powder samples, and thus are not sensitive enough to 
the weak structural/magnetic order to allow a conclusive 
determination on the nature of the structural and AF 
phase transitions near optimal superconductivity. 

In this Letter, we report systematic transport, NMR, 
x-ray and neutron scattering studies of BaFe 2 (Asi_ x P ;E )2 
single crystals focused on determining the P-doping 
evolution of the structural and magnetic phase tran¬ 
sitions near x = 0.3. While our data for x < 0.25 
are consistent with the earlier results obtained from 
powder samples [34], we find that nearly simultaneous 
structural and magnetic transitions in single crystals of 
BaFe 2 (Asi_ ;E P ;E )2 occur at T s T/v > T c = 29 K 
for x = 0.28 and 0.29 (near optimal doping) and dis¬ 
appear suddenly at x > 0.3. While superconductivity 
dramatically suppresses the static AF order and lattice 
orthorhombicity below T c for x = 0.28 and 0.29, the 
collinear static AF order persists in the superconduct- 
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FIG. 2: (a,c,e) Wave vector scans along the [H : 0, 3] direction 
at different temperatures for x — 0.19,0.28, 0.29, and 0.31, 
respectively. Horizontal bars indicate instrumental resolution. 
(b,d,f) Temperature dependence of the magnetic scattering 
at Qaf = (1,0,3) for x — 0.19,0.28, and 0.29, respectively, 
(g) NRSE measurement of temperature dependence of the 
energy width (T is the HWHM of scattering function and zero 
indicates instrumental resolution limited) at Qaf = (1,0,3) 
for x = 0.29. (h) The magnetic order parameters from the 
normal triple-axis measurement on the same sample. 


ing state. Our neutron spin echo and NMR measure¬ 
ments on the x = 0.29 sample reveal that only part of 
the sample is magnetically ordered, suggesting its meso¬ 
scopic coexistence with superconductivity. Therefore, in 
spite of reduced impurity scattering, P-doped BaFe 2 As 2 
has remarkable similarities in the phase diagram to that 
of electron-doped Ba(Fei_ x T c ) 2 As 2 iron pnictides with 
an avoided QCP. 

We have carried out systematic neutron scatter¬ 
ing experiments on BaFe 2 (Asi_ a ,P a ,)2 with x = 
0.19,0.25,0.28,0.29,0.30, and 0.31 [37] using the C5, 
RITA-II, and MIRA triple-axis spectrometers at the 
Canadian Neutron Beam center, Paul Scherrer Institute, 
and Heinz Maier-Leibnitz Zentrum (MLZ), respectively. 
We have also carried out neutron resonance spin echo 
(NRSE) measurements on the x = 0.29 sample using 
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FIG. 3: Temperature evolution of S for (a) x — 0.19 and (b) 
x — 0.28 samples. The solid circles indicate X-ray data where 
clear orthorhombic lattice distortions are seen. The open cir¬ 
cles are data where one can only see peak broadening due to 
orthorhombic lattice distortion. Temperature dependence of 
the [H, 0, 0] scans for (c) x = 0.29 and (d) x — 0.31. The ver¬ 
tical color bar indicates X-ray scattering intensity. The data 
were collected while warming system from base temperature 
to a temperature well above T s . 


TRISP triple-axis spectrometer at MLZ [35]. Finally, 
we have performed high-resolution x-ray diffraction ex¬ 
periments on identical samples at Ames laboratory and 
Advanced Photon Source, Argonne National Laboratory 
[36] . Our single crystals were grown using Ba 2 As 2 /Ba 2 P 3 
self-flux method and the chemical compositions are de¬ 
termined by inductively coupled plasma analysis with 1% 
accuracy m- We define the wave vector Q at (q x ,q y ,qz) 
as (iL, iF, L) = (q x a/ 2tt , q y b / 2tt , q z c /27t) reciprocal lat¬ 
tice units (rlu) using the orthorhombic unit cell suit¬ 
able for the AF ordered phase of iron pnictides, where 
a ~ b ~ 5.6 A and c = 12.9 A. Figure 1(b) shows temper¬ 
ature dependence of the resistivity for x = 0.31 sample, 
confirming the high quality of our single crystals [28]. 

Figure 1(c) summarizes the phase diagram of 
BaFe 2 (Asi_ ;E P ;E )2 as determined from our experiments. 
Similar to previous work on powder samples with x < 
0.25 [34], we find that the structural and AF phase tran¬ 
sitions for single crystals of x = 0.19,0.28, and 0.29 oc¬ 
cur simultaneously within the sensitivity of our measure¬ 
ments (~1 K). On approaching optimal superconductiv¬ 
ity as x 0.3, the structural and magnetic phase transi¬ 
tion temperatures are suppressed to T s « T/v ~ 30 K for 
x = 0.28, 0.29 and then vanish suddenly for x = 0.3, 0.31 
as shown in the inset of Fig. 1(c). Although supercon¬ 
ductivity dramatically suppresses the lattice orthorhom- 
bicity and static AF order in x = 0.28, 0.29, there are still 
remnant static AF order at temperatures well below T c . 
However, we find no evidence of static AF order and lat¬ 
tice orthorhombicity for x = 0.3 and 0.31 at all temper¬ 
atures. Since our NMR measurements on the x = 0.29 
sample suggest that the magnetic order takes place in 


about ^50% of the volume fraction, the coupled T s and 
T/v AF phase in BaFe 2 (Asi_ a ,P a ,)2 becomes a homoge¬ 
neous superconducting phase in the weakly first order 
fashion, separated by a phase with coexisting AF clus¬ 
ters and superconductivity [dashed region in Fig. 1(c)]. 

To establish the phase diagram in Fig. 1(c), we first 
present neutron scattering data aimed at determining 
the Neel temperatures of BaFe 2 (Asi_ a ,P a ,) 2 . Figure 2(a) 
shows scans along the [if, 0,3iL] direction at different 
temperatures for the x — 0.19 sample. The instrumen¬ 
tal resolution limited peak centered at Qaf = (1,0,3) 
disappears at 99 K above T/v [Fig. 2(a)]. Figure 2(b) 
shows the temperature dependence of the scattering at 
Qaf = (1,0,3), which reveals a rather sudden change 
at T/v = 72.5 zb 1 K consistent with the first order na¬ 
ture of the magnetic transition [53] , Figure 2(c) plots 
[if, 0,0] scans through the (1,0,3) Bragg peak showing 
the temperature differences between 28 K (4 K) and 82 
K for the x = 0.28 sample. There is a clear resolution- 
limited peak centered at (1,0,3) at 28 K indicative of 
the static AF order, and the scattering is suppressed but 
not eliminated at 4 K. Figure 2(d) shows the tempera¬ 
ture dependence of the scattering at (1,0,3), revealing 
a continuously increasing magnetic order parameter near 
T/v and a dramatic suppression of the magnetic inten¬ 
sity below T c . Figures 2(e) and 2(f) indicate that the 
magnetic order in the x = 0.29 sample behaves simi¬ 
lar to that of the x = 0.28 crystal without much reduc¬ 
tion in T/v- On increasing the doping levels to x = 0.3 
[36] and 0.31 [Fig. 2(f)], we find no evidence of mag¬ 
netic order above 2 K. Given that the magnetic or¬ 
der parameters near T/v for the x = 0.28, 0.29 samples 
look remarkably like those of the spin cluster phase in 
electron-doped Ba(Fei_ x T c ) 2 As 2 near optimal supercon¬ 
ductivity [22j [23], we have carried out additional neutron 
scattering measurements on the x = 0.29 sample using 
TRISP, which can operate as a normal thermal triple¬ 
axis spectrometer with instrumental energy resolution of 
A E ~ 1 meV and a NRSE triple-axis spectrometer with 
A E « 1 /ieV [35] • Fig. 2(h) shows the triple-axis mode 
data which reproduces the results in Fig. 2(f). How¬ 
ever, identical measurements using NRSE mode reveals 
that the magnetic scattering above 30.7 K is quasielas¬ 
tic and the spins of the system freeze below 30.7 K on a 
time scale of r ^ h/AE « 6.6 x 10 -10 s [23]. This spin 
freezing temperature is almost identical to those of nearly 
optimally electron-doped Ba(Fei_ a ,T c ) 2 As 2 [2IH23] . 

Figure 3 summarizes the key results of our x-ray scat¬ 
tering measurements carried out on identical samples 
as those used for neutron scattering experiments. To 
facilitate quantitative comparison with the results on 
Ba(Fei_ a ,T c ) 2 As 2 , we define the lattice orthorhombic¬ 
ity S = (a — b)/ {a + b) jT9[ [22] . Figure 3(a) shows the 
temperature dependence of 5 for BaFe 2 (Asi_ ;r P ;r )2 with 
x = 0.19, obtained by fitting the two Gaussian peaks in 
longitudinal scans along the (8, 0, 0) nuclear Bragg peak 
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FIG. 4: (a) Temperature dependence of the paramagnetic 
spectral weight for x = 0.25 and 0.29 samples from NMR 
measurements. For x — 0.25, there are no paramagnetic phase 
below 40 K, suggesting a fully magnetic ordered phase. At 
T c of the x = 0.29 sample, there are still 50% paramagnetic 
phase suggesting the presence of magnetic signal outside of 
the radio frequency window of the NMR measurement. The 
spectral weight loss below T c is due to superconductivity. The 
vertical dashed lines mark Tn determined from neutron scat¬ 
tering. (b) The P-doping dependence of the M 2 estimated 
from normalizing the magnetic Bragg intensity to weak nu¬ 
clear peaks assuming 100% magnetically ordered phase. The 
blue solid circles are from [33] • The P-doping levels for dif¬ 
ferent experiments are normalized by their TV’s. The inset 
shows the expanded view of M 2 around optimal doping above 
(solid squares) and below (open squares) T c . (c) The P-doping 
dependence of <5, where the blue diamonds and green squares 
are from Refs. [34] and [28], respectively. For samples near 
optimal superconductivity, the fill and open red circles are S 
above and below T c , respectively. 


[36] . We find that the lattice orthorhombicity S exhibits 
a first order like jump below T s = 72.5 K consistent with 
previous neutron scattering results [34j [36]. We also note 
that the lattice distortion value of S ~ 17 x 10 -4 is similar 
to those of Ba(Fei_ a ,T a ,) 2 As 2 with T s « 70 K ]19l [22]. 

Figure 3(b) shows the temperature dependence of S 
estimated for the x = 0.28 sample. In contrast to the 


x = 0.19 sample, we only find clear evidence of lattice 
orthorhombicity in the temperature region of 26 < T < 
32.5 K [filled circles in Fig. 3(b)] [36] . The open symbols 
represent S estimated from the enlarged half width of sin¬ 
gle peak fits [36]. Although the data suggests a re-entrant 
tetragonal phase and vanishing lattice orthorhombicity at 
low temperature, the presence of weak collinear AF order 
seen by neutron scattering [Figs. 2(c) and 2(d)] indicates 
that the AF ordered parts of the sample should still have 
orthorhombic lattice distortion mm- Figure 3(c) and 
3(d) shows temperature dependence of the longitudinal 
scans along the [i7, 0,0] direction for the x = 0.29 and 
0.31 samples, respectively. While the lattice distortion 
in the x = 0.29 sample behaves similarly as that of the 
x = 0.28 crystal, there are no observable lattice distor¬ 
tions in the probed temperature range for the x = 0.31 
sample. 

To further test the nature of the magnetic ordered state 
in BaFe 2 (Asi_ a ,P a ,) 2 , we have carried out 31 P NMR mea¬ 
surements under a 8-T c-axis aligned magnetic field |36j. 
Figure 4(a) shows the temperature dependence of the in¬ 
tegrated spectral weight of the paramagnetic signal, nor¬ 
malized by the Boltzmann factor, for single crystals with 
x = 0.25 and 0.29. For x = 0.25, the paramagnetic spec¬ 
tral weight starts to drop below 60 K and reaches zero 
at 40 K, suggesting a fully ordered magnetic state below 
40 K. For x = 0.29, the paramagnetic to AF transition 
becomes much broader, and the magnetic ordered phase 
is estimated to be about 50% at T c = 28.5 K. Upon 
further cooling, the paramagnetic spectral weight drops 
dramatically below T c because of radio frequency screen¬ 
ing. We find that the lost NMR spectral weight above T c 
is not recovered at other frequencies, suggesting that the 
magnetic ordered phase does not take full volume of the 
sample similar to the spin-glass state of Ba(Fei_ ;E T ;E ) 2 As 2 

Figure 4(b) shows the P-doping dependence of the or¬ 
dered moment squared M 2 in BaFe 2 (Asi_ ;E P ;E )2 includ¬ 
ing data from Ref. [34] . While M 2 gradually decreases 
with increasing x for x < 0.25, it saturates to M 2 ~ 
0.0025 /.l 2 b at temperatures just above T c for x = 0.28 
and 0.29 before vanishing abruptly for x > 0.30. The 
inset in Fig. 4(b) shows the P-doping dependence of the 
M 2 above and below T c near optimal superconductivity. 
While superconductivity dramatically suppresses M 2 , it 
does not eliminate the ordered moment. Figure 4(c) 
shows the P-doping dependence of S in BaFe 2 (Asi_ a ,P a ,)2 
below and above T c . Consistent with the P-doping de¬ 
pendence of M 2 [Fig. 4(b)] and Tn [Fig. 1(c)], we find 
that S above T c approaches to ^ 3 x 10 -4 near optimal 
superconductivity before vanishing at x > 0.3. 

Summarizing the results in Figs. 2-4, we present the 
refined phase diagram of BaFe 2 (Asi_ :E P :E )2 in Fig. 1(c). 
While the present phase diagram is mostly consistent 
with the earlier transport and neutron scattering work on 
the system at low P-doping levels ESDI EH, we have dis- 
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covered that the magnetic and structural transitions still 
occur simultaneously above T c for x approaching optimal 
superconductivity, and both order parameters vanish at 
optimal superconductivity with x = 0.3. Since our NMR 
and TRISP measurements for samples near optimal su¬ 
perconductivity suggests spin-glass-like behavior, we con¬ 
clude that the static AF order in BaFe 2 (Asi_ a ,P x )2 disap¬ 
pears in the weakly first order fashion near optimal super¬ 
conductivity. Therefore, AF order in phosphorus-doped 
iron pnictides coexists and competes superconductivity 
near optimal superconductivity, much like the electron- 
doped iron pnictides with an avoided QCP. From the 
phase diagrams of hole-doped Bai_ IC A c Fe 2 As 2 [SHED, it 
appears that a QCP may be avoided there as well. 
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Supplementary information 

Section A: Details of the neutron and X-ray 
scattering, and NMR experiments 

Neutron scattering experiments: 

We have aligned the x = 0.19 samples in the [ H , A, 3 H] 
scattering plane and the x = 0.25, 0.28, 0.29, 0.30, 0.31 
samples in the [. H , 0, L\ zone.For neutron scattering mea¬ 
surements of the x =0.19 compound at C5 spectrome¬ 
ter, we used a vertically focused PG(002) monochrome- 
tor and a flat PG(002) analyzer with a fixed final energy 
Ef =14.56 meV. We useda PG filter after the sample to 
eliminate the higher order neutrons.At RITA-II, we use 
a PG filter before the sample and a cold Be-filter after 
the sample with the final neutron energy fixed at 4.6 
meV. For MIRA measurements, the final energy was set 
to Ef = 4.06 meV and a Be-filter was additionally used 
as a filter. 

In addition to usual neutron diffraction experiments, 
we have also carried out measurements on TRISP at 
MLZ, Germany. The experimental set these measure¬ 
ments are described in detail in Ref. [17] of the main 
text. 

X-ray scattering experiments: 

The high resolution X-ray diffraction of the x =0.28 


sample was performed using a four circle diffractome¬ 
ter and Cu K a i X-ray radiation from a rotating anode 
X-ray source at Ames Lab. We have used beamline 6- 
ID-D at the Advanced Photon Source at Argonne Na¬ 
tional Laboratory with 100.2 keV incident photon beam 
for measurements of the x = 0.19, 0.25, 0.29, 0.30, 0.31 
compounds.The NMR measurements were performed by 
the Spin-echo technique, and the paramagnetic spectral 
weights were obtained by integrating the spectral in¬ 
tensity at the resonance frequency of the paramagnetic 
phase. 

Section Badditional transport, X-ray and neu¬ 
tron scattering data: 

We have carried careful transport measurement us¬ 
ing 4 probe method in a physical property measurement 
system. Our systematic measurements of the resistiv¬ 
ity for different doping concentrations are shown in Fig. 
SI. Typical raw data for X-ray scattering experiments is 
shown in Fig. S2 for the x =0.19 and 0.28. The presence 
of two peaks along the [iL, 0, 0] direction is a direct in¬ 
dication of lattice orthorhombicity. Figure S3 shows the 
temperature dependence of the magnetic order param¬ 
eter for x = 0.19, 0.25, 0.28, 0.29, 0.30, 0.31 samples. 
Figure S4 shows the raw 31 P NMR spectra for the x = 
0.25 and 0.29 samples at different temperatures. 
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FIG. SI: Temperature dependence of the normalized resistiv¬ 
ity for BaFe 2 (Asi_ x P x ) 2 . The measurements were conducted 
by a standard four-terminal method using a Quantum Design 
Physical Property Measurement System. The inset shows the 
expanded data below 50K. 



FIG. S2: Temperature evolution of the orthorhombic Bragg 
peaks for BaFe 2 (Asi_ cc P cc )2 with x =0.19 and 0.28. The data 
were collected while warming the system from the base tem¬ 
perature. 
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FIG. S3: Temperature dependence of the magnetic order pa¬ 
rameter, where the intensity of the magnetic scattering is ob¬ 
tained by subtracting the data well above TN as background 
and normalized to weak nuclear Bragg peaks.The intensity of 
the x — 0.25 compound was expanded by 4 and x = 0.28, 0.29, 
0.30, 0.31 compounds was expanded by 20. Arrows indicate 
the T c of different samples.The slight intensity differences for 
the x — 0.29 and 0.30 samples are within the error of our 
measurements. The observed magnetic peaks are resolution 
limited, giving an estimated spin-spin correlation length of 
300 A. 
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FIG. S4: The 31 P spectra at different temperatures for x — 
0.25 and x — 0.29 samples. Tn and T c marks the Neel tem¬ 
perature and the superconducting transition temperature, re¬ 
spectively. The horizontal axes show the relative frequency 
from the fixed frequency /o. 
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